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TECHNICAL MEMORANDUM 


A METHOD FOR INCORPORATING CHANGING STRUCTURAL CHARACTERISTICS 
DUE TO PROPELLANT MASS USAGE IN A LAUNCH VEHICLE ASCENT SIMULATION 

1. INTRODUCTION 


Launch vehicles consume large quantities of propellants and oxidizers quickly, causing the mass 
properties and structural dynamics of the vehicle to change dramatically. The current method of handling 
this change for structural load assessments is to simply have a large collection of structural models 
representing various propellant fill levels. This creates a large database of models and compounds the 
problem of having easily derivable reduced models. The offline processing required to generate numer- 
ous models to populate the database can be large and would have to be repeated for model changes. 


2. PROPELLANT MASS METHODOLOGY 


Section 2 presents a method for handling these mass changes in a more efficient manner. 

The method essentially allows for the subtraction of propellant mass as the propellant is used 
in the simulation. This subtraction is done in the modal domain of a launch vehicle generalized model. 
The additional computation involved is primarily in terms of constructing the used propellant mass 
matrix from an initial propellant model and further matrix multiplications and subtractions. There 
is an additional eigenvalue solution required to uncouple the new equations of motion; however, 
this is a much simpler calculation since it starts from a system that is already substantially uncoupled. 

Assume a vehicle model empty of propellants: 


KIRRMRRRJ. 

(1) 

where 


| A/1 = mass matrix of empty vehicle 
| F 1 = stiffness matrix of empty vehicle 

|x y } = acceleration of empty vehicle degrees of freedom (dofs) 

jx v j = displacement of empty vehicle dofs 

{F v } = applied forces acting on the empty vehicle, 


reduced modally: 


R} = MR) 

(2a) 

R} = MR} 

(2b) 

iR«d + K]R} = MR} • 

(3) 


where 


= eigenvectors or modeshapes of unconstrained empty vehicle 
{ <7 v } = generalized displacement of empty vehicle generalized dofs 
j<r/ v ,] = generalized acceleration of empty vehicle generalized dofs 
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[/] = identity matrix resulting from mass normalized modeshapes 



= eigenvalues of unconstrained empty vehicle. 


Now, assume a propellant load on a subset of the original vehicle dofs: 


[M p ]jx v/5 ] , 


( 4 ) 


where 


[M p \ = mass matrix of propellant model 

{i vP }= accelerations of propellant model dofs selected from |i v j. 


Using the same modal transformation, and selecting the appropriate rows for the propellant 
dofs, [<£> J, the propellant mass can be reduced: 


where 


[® v p\ [ M p\ f^vl R} 



( 5 ) 


[**] = propellant model selected rows from [0 y j 
= propellant model mass matrix converted to 


empty vehicle generalized dofs. 


M 


is not necessarily diagonal because the modes are no 


model can now be coupled to the vehicle: 


longer orthogonal; however, this propellant 





(6) 


This model can then be reduced modally once again. This time the reduction is at a much reduced 
cost because of the reduced number of dofs and the uncoupled nature of the nonpropellant portions 
of the equation: 


(9.1 = [4*2] {9,2} 

(7a) 

{«>■} = Nfe} 

(7b) 

1T/A2} + [^2] {912} = [^2] [k] k} * 

(8) 


where 


[<Z> 2 ] = eigenvectors of modally modified unconstrained empty vehicle 

R 2 } = generalized displacement of modally modified empty vehicle generalized dofs 
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fe} 

[/] 


= generalized acceleration of modally modified empty vehicle generalized dofs 
= identity matrix resulting from mass normalized modeshapes 

= eigenvalues of modally modified unconstrained empty vehicle. 


This would once again make the modes orthogonal and uncouple the equations of motion for easy 
solution. 


If all modes are kept during this process, it is easy to see that the solutions would be the same. 
However, if modes are truncated, a good deal of information could be lost. This is because the empty 
vehicle structure would have higher frequency modes that could be truncated. In reality, some of these 
modes would be lowered in frequency by the presence of the propellant mass and should be retained. 


This truncation problem can be compensated for in the following way. The first vehicle model 
is constructed for a fully loaded vehicle with propellant model. This would assure that all the lower fre- 
quencies of interest are retained. The subset propellant model is then actually a model of the propellant 
used or removed. Equation (6) would then become 



M 


,]){«,} + K1w=mw- 


(9) 


Note the change in sign in the first term. This is to reflect that mass is being removed from the vehicle. 


Equation (9) can then be reduced further as in equation (8). This uncoupled set of equations 
of motion could then be solved to provide the response of the complete system. Since the equations 
are uncoupled, they could be solved using a closed-form solution of a second-order dynamic equation 
with initial conditions. This solution still assumes the mass and the dynamics of the system remain con- 
stant during the time of solution. However, there is now an approach that would allow a much greater 
frequency of updating the mass characteristics of the system so the approximation has less of an effect. 


Load transformation matrices (LTMs) used for data recovery are most easily defined in terms 
of the initial model governed by equations (1) and (3). For the displacement method of data recovery, 
generating the correct response is simply a matter of including the additional transformation of 
equation (7): 


{Resp} = [LTM][0 v ]{q v } = [LTM] [<Z> V ] {<Z>2 }{q v2 } , 

where 


( 10 ) 


{Resp}= desired system response 
[LTM] = displacement method LTM. 
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For acceleration method LTMs, some additional work is required because of the inclusion of an 
applied force term. Equation (9) can be rewritten by moving the propellant mass term to the right-hand 
side and treating it as an applied load: 




( 11 ) 


In this case, the acceleration term simply includes the equation (7) transformation. The applied force 
term must include the propellant inertial load term: 


{Resp} = [LTMA][0 v ][02\{q v2 \ + [LTMD] [<Z> v ] [02] {q v2 } 

+[LTMF] ([«>„]{?„} + [M p ][4>2]{^ 2 }) , (12) 

where 


{Resp} = desired system response 
[LTMA] = acceleration term of acceleration method LTM 
[LTMD] = displacement term of acceleration method LTM 
[LTMF] = applied force term of acceleration method LTM. 
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3. EARLY PROTOTYPE 


In 2001, the Space Launch Initiative program requested the Marshall Space Flight Center engi- 
neering groups to study what would be required to perform integrated system-level technical assess- 
ments of candidate launch systems. From this, the Vehicle Integrated Performance Analysis team was 
formed as a prototype answer and culminated in an assessment of the Saturn V ascent vehicle. This 
vehicle was chosen to validate the process because of the availability of data and the lack of proprietary 
or political issues. 

Figure 1 shows a comparison of first-stage bending moments for a given flight loading condition 
with varied propellant mass. The solid lines indicate results with the propellant masses modeled seper- 
ately for each propellant fill level. The point marker curves indicate results using the full propellant mass 
modified by the mass reduction method described in this Technical Memorandum. Further investigation 
needs to take place regarding modal truncation effects; however, the results are excellent. The full 
propellant mass vehicle model was reduced carrying modes up to ~25 Hz. 
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Figure 1. Mass model results comparison; VAC-02 bending moment— max Q-alpha loads. 
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4. CONCLUSION 


Launch vehicles consume large quantities of propellants and oxidizers quickly. A method 
has been presented for handling these mass changes more efficiently for a structural load analysis 
simulation. The method allows for the subtraction of propellant mass as the propellant is used. 

The subtraction is done in the modal domain of a launch vehicle generalized model. The method 
has been successfully tested in a simulation of Saturn V loads. Results from the use of the method 
have been compared to results from separate structural models for several propellant levels. This 
comparison shows excellent agreement. Effects of modal truncation should be assessed when using 
the method. The method could also be further developed to encompass more complicated propellant 
models, including slosh dynamics. 


7 


REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operation and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 

1. AGENCY USE ONLY (Leave Blank) 

2. REPORT DATE 

November 2004 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 



5. FUNDING NUMBERS 

A Method for Incorporating Changing Structural Characteristics 
Due to Propellant Mass Usage in a Launch Vehicle Ascent Simulation 


6. AUTHORS 

D.S. McGhee 


7. PERFORMING ORGANIZATION NAMES(S) AND ADDRESS(ES) 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, AL 35812 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

M-1124 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-000 1 

10. SPONSORING/MONITORING 
AGENCY REPO NUMBER 

NAS A/TM - 2004-2 1 3549 

11. SUPPLEMENTARY NOTES 




Prepared by the Engineering Directorate 



12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 

Unclassified-Unlimited 
Subject Category 15 

Availability: NASACASI 301-621-0390 



13. ABSTRACT (Maximum 200 words) 

Launch vehicles consume large quantities of propellant quickly, causing the mass properties and structural dynam- 
ics of the vehicle to change dramatically. Currently, structural load assessments account for this change with a large 
collection of structural models representing various propellant fill levels. This creates a large database of models 
complicating the delivery of reduced models and requiring extensive work for model changes. Presented here is 
a method to account for these mass changes in a more effiicient manner. The method allows for the subtraction 
of propellant mass as the propellant is used in the simulation. This subtraction is done in the modal domain of the 
vehicle generalized model. Additional computation required is primarily for constructing the used propellant mass 
matrix from an initial propellant model and further matrix multiplications and subtractions. An additional eigen- 
value solution is required to uncouple the new equations of motion; however, this is a much simplier calculation 
starting from a system that is already substantially uncoupled. The method was successfully tested in a simulation 
of Saturn V loads. Results from the method are compared to results from separate structural models for several 
propellant levels, showing excellent agreement. Further development to encompass more complicated propellant 
models, including slosh dynamics, is possible. 

14. SUBJECT TERMS 

structural loads, simulation, mass properties, propellant usage, 

15. NUMBER OF PAGES 
12 

launch vehicle, modal reduction, model reduction, structural dynamics 

16. PRICE CODE 

17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 

20. LIMITATION OF ABSTRACT 

Unlimited 

NSN 7540-01-280-5500 



Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 


298-102 


